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ABSTRACT 

Context. We present a theoretical investigation of dielectronic recombination (DR) of Ar-like ions that sheds new light on the behavior 
of the rate coefficient at low-temperatures where these ions form in photoionized plasmas. 

Aims. We provide results for the total and partial Maxwellian-averaged DR rate coefficients from the initial ground level of Kn- 
Znxm ions. It is expected that these new results will advance the accuracy of the ionization balance for Ar-like M-shell ions and pave 
the way towards a detailed modeling of astrophysically relevant X-ray absorption features. 

Methods. We utilize the AUTOSTRUCTURE computer code to obtain the accurate core-excitation thresholds in target ions and carry 
out multiconfiguration Breit-Pauli (MCBP) calculations of the DR cross section in the independent-processes, isolated-resonance, 
distorted-wave (IPIRDW) approximation. 

Results. Our results mediate the complete absence of direct DR calculations for certain Ar-like ions and question the reliability of 
the existing empirical rate formulas, often inferred from renormalized data within this isoelectronic sequence. 

Key words. Atomic data - Atomic processes 



1. Introduction 

Atomic structure and dynamic behavior of highly-charged ions 
is one of the key ingredients presently required for both lab- 
oratory plas ma diagnostics and interpreta t ion of astrophysical 
pheno mena dFawcettl Il99ll iLiedhall l2000t iKallman & Palmeril 
I2007I) . The most common diagnostic technique to probe for elec- 
tron temperature of laser-produced, fusion, or astrophysical plas- 
mas involves spectroscopical observations of intensity ratios of 
EUV or X-ray emission lines coming from consecutive ioniza- 
tion stages of a single plasma component. The accuracy of de- 
rived plasma parameters is strongly affected by uncertainties in 
chemical abundances, often inherited through the use of unre- 
liable collis ional ionization and/or dielectronic recom bination 
(DR) rates dSavin & Laming! |2002|) . ISeon et all d2003l) investi- 
gated the effect of the uncertainties in DR rates on an isoelec- 
tronic line ratio in Ti and Cr plasmas and found substantial dif- 
ferences in the fractional abundances obtained for Ti v and Cr vn 
ions (a shift of the curves to lower temperatures) as compared to 
those in ferred using the recom mended recombination rate coef- 
ficients dMazzotta et al.lll998l) . 

In recent years an enormous amount of progress, both the- 
oretically and experimentally, has been made in improving the 
DR rate database along isoelectronic series within the first and 
second ro ws, and the third row up th rough Mg-like ions (see the 
review bv lKallman & Palmeril (I2007I) . and references therein). In 
most cases it is found that the newly determined DR rates are sig- 
nificantly larger than their earlier recommended estimates, hav- 
ing profound consequences on the ionization balance and ther- 
mal equilibrium in both photoionized and collisionally ionized 
plasmas, from the solar corona to Ac t ive Galactic Nuclei (AGN) 
(see, for exam pl e. |Chakravortv et al.l J2008. 2009l): lB~rvans et ail 
(l2006ll2009allbl) : lDere et al.l(l2009h l In this work, we present im- 
proved theoretical predictions of the DR rates of Ar-like ions. Of 
these ions, those of the iron peak elements are some of the more 



abundant in cosmic plasmas. In much of the remainder of this in- 
troduction we will provide some of the motivation in improving 
their atomic database, in particular their DR rate coefficients. 

The strong contribution of M-shell Fe ions to the unresolved 
transition array of inner-shell absorption lines in ~ 15- 17 A X- 
ray spectra of several AGN observed with XMM-Newto n and 
Chandra, was initially not well understood dNetzerl 120041) . The 
fact that AGN photoionization models initially overpredicted the 
average ionization stage of iron was attributed in part to an un- 
derestimat e in the low-tempera ture DR rate coefficients for M- 
shell iron dKraemer et al J 120041) . and constituted the main mo- 
tiv ation behind the benchmark calculations recently performed 
bv lBadnelll d2006h as well as exper imental and theoretical re- 



sults presented in 



Lukic et al 



_J 2007b. The present work further 
extends the calculations of Badnelll d2006l) for Fe ix ions by aug- 
menting the configuration interaction (CI) with some of the most 
important Aw c = ionic core excitations dAggarwal et al.|[2006h 
IZeng et alJl2006l) . 

In the f ramew ork of testing nucleosynthesis models, 
lEllison et al ] d2001l) identified Co as a rewarding element to 
study galactic and stellar formation histories th rough the ob- 
served abundance trends ddel Peloso et al.ll2005l) . For example, 
the agreement of the modeled time-dependent ejecta composi- 
tions and velocities with observed Type la supernovae spectra 
requires a substantial initial p resence of Ni, Co, and Fe i n 
the outer layers of the ejecta dHillebrandt & Niemeverll2000l) . 
In the present work, we study the photorecombination of Co x 
as an example of a heavier iron peak element for which the 
stellar photospheric e lemental abundances are less well known 
dAdelman et al.ll2000t) . 

Nickel is one of the most imp ortant heavy impurit ies in toka- 
maks and early attempts by the ITFR Groupl ((1980) to model 
the fractional abundance of its charge states suffered from de- 
ficient DR data. Recent simulations of the observed plasma 
emission from magnetic confinement fusion devices, namely the 
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JET t okamak in Abingdon and RFX in Padova dMattioli et alJ 
12004 . also lacked accurate recombination rate coefficients for 
the Ni xi ion. However, reliable electron impact ionization data 
of the remaining ions in the nickel iso nuclear sequence have 



as 



be en provided by Pindzola et al.l d!991l) and were widely used 
bv lMattioli etalT(l2004l) in simulations of Ni emission line spec- 
tra. The past few years have marked a renewed theoretic al inter- 
est dVerma et al.ll2007l:lAggarwal & Keenanl2007ll2008l) in elec- 
tron excitation data for argon-like nickel, initiated by its identi- 
ficatio n in numerous astrophysical plasmas (consult rVerma et all 
(120071) for an exhaustive up-to-date bibliography). The most 
recent studie s of th e Intra-Cluster Medium, as discussed by 
IWerner et al.l d2008l) . put constraints on supernova models by 
using Ni/Fe abundance patterns in the ejecta of type la super- 
novae. In addition, from the Mass Time-of-Flight Spectrometer 
data accumulated d uring the first decade of SOHO's operation, 
iKarrer et al.l d2007l) inferred charge-state distributions, isotopic 
composition, and the elemental Ni /Fe ratio of the solar wind, 
and confirmed that both nickel and iron become enriched in the 
solar corona. 

The si gnificance of rel i able a tomic data has been demon- 
strated by IChurazov et al. I J200l through the modeling and 
interpretation of the 5-9 keV spectrum from the multi- 
temperature cor e of th e Perseus galax y cluster. In their study, 
IChurazov et"aTI d2004l) used APEC Ismith et all 1200 lb and 
MEKAL dMewe et alii 19851) models, both having the redshift 
of major line energies, the heavy element abundances, and the 
plasma temperatures as free parameters. It has been shown 
that the MEKAL model yields the best description of the 
spectra only when nickel is overabundant relative to iron by 
a factor of ~ 2 compared to solar. However, this enhance- 
ment is not required by the simulations of the APEC code 
(http://cxc.harvard.edu/atomdb) that uses updated atomic data. 

The present computational study is part of a n ongoing in- 
vestigation of PR processes in argon-like ions dNikolic et alJ 
120071 120091 l2010h . and deals with A« c = 0, 1 ionic core 
excitations and associated dielectronic resonances that domi- 
nate electron-ion recombination in photoionized plasmas. The 
theoretical foundation and computational method we use for 
the DR calculations are found elsewhere dBadnelletalJl200l 
and here we only outline the essence. With the use of th e 
open-source AUTOSTRUCTURE code dBadnelllll986l Il997l) . 
we carry out MCBP computations of energy levels and decay 
rates in an intermediate coupling scheme for Ar-like K n, Ca in, 
Sciv, Tiv, Vvi, Crvn, Mnvm, Feix, Cox, Nixi, Cuxn, and 
Znxm ions. In order to account efficiently for all DR contri- 
butions coming from numerous Rydberg series of resonances 
and offer them to the plasma modeling community in conve- 
nient level-resolved format, we further enforce the independent- 
processes, is olated-resonance, d istorted-wave (IPIRDW) ap- 
proximation dPindzola et al.lll992l) . 

In the remainder of this paper, we will proceed as follows. 
Section [2] discusses the process of electron-ion recombination 
and transparently outlines the main relations and equations aris- 
ing within the adopted methodology. A comparative overview of 
existing atomic structure in argon-like ions is provided through- 
out Section [3] and an analysis of the results is presented in 
Section[4] 



2. Elementary processes of relevance 

The contribution of the photorecombination process of an ion- 
ization state q+ from a single partial wave (J, n) can be described 
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Here c is a continuum state consisting of an initial electron in- 
cident upon the target ion I (9)+ that is either directly captured, 
via radiative recombination (RR), to a I (f ^ 1,+ * bound state, b, 
or captured into an autoionizing doubly-excited state, 

d, that undergoes subsequent radiative decay to the same final 
bound state, b, completing the DR process. The present work 
investigates {3s, 3p] — > {3d, At'} inelastic excitations from the 
ground state that give rise to Rydberg series of autoioniz- 
ing states formed by the capture of the scattered electron. Each 
Rydberg series of resonances will converge to a corresponding 
threshold, given in Tables [3] |4]|5]|6] and [7] 

Within the adopted IPIRDW framework, our calculations 
rely on lowest-order perturbation theory to compute Lorentzian 
resonance profiles as a function of the electron's center-of-mass 
(cm.) energy e. The total DR cross section is then given as 



c n ^ ^ ^ 
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Here ej = 2(e - e</)/r</ is the reduced cm. energy at which the 
(partial) integrated resonant strength is 



S^JMbRy 2 ] = 2.674- 10" 
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with gd and g, „ as statistical weights of the resonant state d in 
the recombined ion and the ground state of the target ion, respec- 
tively. The summations over c' and s in Eq. (f3]) cover all states 
that are attainable from resonant state d either by radiative decay 
or by autoionization, with corresponding rates given in inverse 
seconds. Hence, the sum over s includes not only bound states 
that are below the first ionization limit of the recombined ion, 
Ei?t but also may consider a radiative cascade through other au- 
toionizing states resulting in the total radiative rate, A 1 ,. The sum 
over c' accounts for resonant scattering (excitation) and amounts 
to the total autoionization rate, A a d . In Eqs. <f2j»-([3j we implicitly 

assume that both the resonance position, ed - Ed - EJ, and the 
total width, Fd = h(A a d + A r d ), are in Rydberg units. In addition, 
the summation over d spans all of the N J -" doubly excited states 
of given parity n and total angular momentum / that are formed 
through Eq. dT) for 3 < n < 1000 and < £ < 10, wherein the 
index b considers all accessible bound states. 

The DR/RR rate coefficient (in units of cm 3 s ) for \ q+ 
ions in a plasma with a Maxwellian electron energy distribution 
/mb(c, T) is given by 



Jo 



DR/RR, 



(4) 



and conveniently modeled using physically-motivated fitting for- 
mulae 



a DR (T) = T- yl Y a exp(-Ei/T) and 



(5) 
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a RR (T) 



(6) 



where the more general form p ± = l+B + C exp(-72/r) is used 
here. The fitted DR rate coefficients given by Eq. (0 are rela- 
tively smooth and broad curves with several local maxima, each 
at a particular temperature T,- = 2E//3, with £, given in Table [8] 
On the other hand, the RR rate coefficients given by Eq. (O de- 
crease monotonically with temperature; the fitting coefficients 
A, B, C, and To, 1,2 are given in Table [9] These nonlinear least- 
squares fits are accurate to better than 1% over (10 1 - 10 7 )q 2 K 
with the correct asymptotic forms outside of this temperature 
range. 



3. Atomic structure 

Before performing the electron-ion scattering calculation in 
Eq. ([TJ, it is necessary to obtain an accurate description of the 
Ar-like target-ion states. Our approach for the entire Ar-like iso- 
electronic seq uence is essentially the same as what was done 
earlier for Ti v dNikolic et al.l2 009). Since the dominant DR con- 
tribution is due to the e~ + 3s 2 3p 6 — > 3s 2 3p 5 3d n£ capture, we 
are particularly concerned with obtaining highly-accurate wave- 
functions for both the 3s 2 3p 6 ground state and the 3s 2 3p 5 3d 
excited state. However, we also need to consider the 3s3p 6 3d, 
3s 2 3p 5 AC , and 3s3p 6 4i' target states that also contribute to the 
DR rate coefficient via Eq. ([]]). Single and double promotions 
out of these configurations are also included in our total config- 
uration expansion for each target state. 

With this configuration basis, the Is, 2s, 2p, 3s, 3p, and 3d 
Slater-type orbitals were determined by varying the radial scal- 
ing parameters, A„c, to minimize the equally-weighted-sum of 
MCBP eigenener gies of the seventeen lowest states, using the 
model potential of lBurgess et alj dl989) with mass-velocity and 
Darwin corrections. The values of the radial scaling parameters 
for closed-core orbitals {Is, 2s, and 2p) remained throughout the 
optimization in near proximity to their default value of 1 .0. We 
then applied small adjustments to the 31 scaling parameters in 
order t o reproduce the N IST experimental 3p — > 3d radiative 
data o f IShirai et alj J2000h . The resultant A„c values are listed in 
Table [1] and the radiative data are listed in Table [2] As is seen 
in Table [2 our target description is such that the present oscil- 
lator strengths are in excellent agreement with the experimental 
values. 

We have also confirmed that the computed eigenenergies 
were stationary with respect to the small variations in A n c val- 
ues that were used to fine-tune the 3p — > 3d oscillator strengths. 
Indeed, as is seen in Tables [3] |4j |5J |6l and [7] for the lowest- 
lying Ar-like states, our computed energies are in excellent 
agreement with the exper imental values given in the NIST 
Atom i c Spectra Database dSugar & Corliss! 119851; IShirai et aT] 
l2000t iPettersen et alJ 120071) . In the c ase of higher ionization 
stages , the atomic data produced by llrimia & Froese-Fischeri 
using the multi-configurational Hartree-Fock (MCHF) 
m ethod, or by a non-re l ativist ic single-configuration approach 
of iGhosh & Mukherieel dl997l) . are up to 7.2% lower than the 
NI ST values. The bas ic CIV3 atomic data for Feix produced 
by IVerma et al.l (|2Q06|) are noticea b ly clo ser to the present re- 
sults than those of lAggarwal et al.l (120061) using the fully r ela- 
tivistic multi-configurational Flexible Atomic Code (FAC) of 
(12004 . 



4. DR Results 

To treat the DR processes occurring in Eq. ([TJ, all possible con- 
tinuum and resonance wavefunctions are constructed by cou- 
pling an appropriate distorted-wave free et or bound nt orbital 
to each target configuration wavefunction, as obtained in the pre- 
vious section. The energies, radiative rates, and Auger rates of 
each resonance are then computed and used in Eqs. (0 and (0| 
to produce Maxwellian rate coefficients that are shown in Figs. [7] 
[8] and [9] In all the figures, the bars indicating the collision- 
ally ionized zone were obtained assuming the conditions of 
coronal equilibrium with electron-collisional plasma tempera- 
tures for which the fractional abundance of the target ion i n 
question surpasses 10% of its peak value dBrvans et al.ll2008h . 
Similar approximate temperature indicators for gas in photoion- 
ization equilibrium have been computed using Cloudy (v08.00; 
iFerland et al.l (ll998l) V Although these results did not include the 
DR rate coefficient data reported here, the indicated temperature 
range should still be reasonably accurate. 

For K 11 and Ca m, the only o the r data available in Fig. [TJ are 
the results of lMewe et~aLl(ll980h andlMazzotta et al.l dl998l) . The 
rate coefficients of Mazzotta et al. J 1998b are inferred by using 
the general formula of lBurgessI (119651)" that is first scaled down to 
match the empirical DR rate of lHahnl (Il989l) for Fe ix and then 
this same scaling is applied to other ions in the isoelectronic 
sequence. Surprisingly, this procedure gives results that are in 
fair agreement with our present calculations, but, as is seen for 
higher ionization stages, this method fa ils to account corre ctly 
for low-temperature DR. The results of iMewe et a i] (ll980l) . on 
the other hand, were obtained using a parameter ized empiri- 
cal for mula based upon the reno rmalized results of iJacobs et al.l 
d 1977b and lAnsari etaD dl970h for 3p -> 3d, 3p -> 4s, and 
3p — > Ad core-excitations of Feix; that procedure gives erro- 
neous results for the position and peak of the rate coefficient. 

The results shown in Fig. [2] have been ana l yzed m ore fully 
in our earlier studi es of Sc iv dNikolic et alJ 1201 Oh and Ti v 
dNikolic et al.l 12009). but the following points should be made. 
First, there is experimental dat a available for both of these 
10ns dSchippers et alJl998Ll2002h . and these measurements were 
useful for quantifying the positions of low-energy resonances. 
However, the experiments for these low-charged Ar-like ions 
were subject to motional Stark effect that reionized the higher-n 
recombined bound states and led to a reduction in the DR cross 
section. Thus, the measured rate coefficient is unphysically too 
low at higher temperatures and is inappropriate for use in plasma 
modeling. Also, the experimental rate coefficient also contains 
the contribution due to RR, which is why the Sc iv experimental 
rate coefficient exceeds our DR result at lower temperatures and 
merges with our RR result. 

Another new feature seen in Fig.[2]that was absent in Fig. UJ 
is that, as the ionization stage is increased, the lowest-lying 
members of the Rydberg series given in Eq. ([TJ, namely the 
3s 2 3p 5 3d 2 and 3s 2 3p 5 3dAs resonances, approach zero contin- 
uum energy and give a large contribution to the rate coefficient at 
the lower temperatures found in photoion ized plasmas. This con - 
tribution is not included in the results of lMazzotta et alj dl998l): 
those data are determined from the Burgess formula dBurgessI 
dl965h ). which only considers high-t emperature DR. We note 
that the data of lMazzotta et al.l (fl998: : ) also differs significantly 
from the present DR rate coefficient in the collisionally ion- 
ized zone. The other available data - the e mpiri cal results of 
IMewe et al.l d!980l) . lLandini & Fo"ssil d!991l) . and iHahnl d!991l) 
- differ dramatically from the present MCBP results. 
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For the higher ionization stages of V vi, Cr vn, and Mn vm 
shown in Figs. |3]and|4] the present rate coefficient again has a 
significant low-temperature DR contribution that is absent in all 
previously available data. For Fe ix, there exist several other re- 
sults, but most of these are again empirical results that fail to 
account for low-temperature DR contributions in the photoion- 
ized plasma region and also differ significantly from our present 
DR results in the collisionally-i onized p l asma region. We also 
show the earlier MCBP results o f lBadnellld2006l) . Those calcula- 
tions were quite similar to the present ones, except that a smaller 
configuration basis was used and thus a slightly larger rate co- 
efficient was obtained. We have established that by augmenting 
those earlier calculations with the inclusion of additional cor- 
relation configurations, as we include here, the two results are 
brought into agreement. Both MCBP results are somewhat lower 
than the expe rimental rate coefficie nt, as has been discussed 
more fully by Sc hmidt et al.l (120081) . We note that the experi- 
ment for Fe ix was not influenced by external fields, unlike the 
experiments for Sc iv and Ti v. Our DR rate coefficients for Co x 
through Zn xm are shown in Figs.[5]and[6] where it is seen again 
that all previous data do not account for the low-temperature 
contributions that dominate in the photoionized plasma region. 

In Figs.|7][8j and|9] we show the DR contributions from each 
resonance series. As anticipated due to its large core oscillator 
strength, the 3s 2 3p 5 3dn( resonances are dominant. The lowest 
3s 2 3p 5 3d 2 resonances are first seen to be indistinct from the rest 
of the series for Kn and Cam, then begin to show a separate 
feature at low temperatures for Sc iv, and completely dominate 
at low temperatures for Ti v. As the nuclear charge is increased 
from V vi through Zn xm, additional low-temperature features 
are seen to oscillate as the n — 4 and n — 5 resonances move 
from above to below threshold. Also seen as the ionization stage 
increases is the appearance of the 3s 2 3p 4 3d 2 n€ (n = 3,4) res- 
onances near threshold that dominate the low-temperature rate 
coefficient. 

Lastly, we list the DR fitting coefficients for each ion, as 
described by Eq. (01, in Table [8] The RR fitting coefficients of 
Eq. (O are listed in Table [9] 

5. Summary 

It is clear that all previously-available DR data for Ar-like ions 
are inadequate at lower temperatures since the contributions 
from low-lying resonances are not accounted for. Thus, those 
earlier data are inappropriate for use in modeling photoionized 
plasmas. And even at higher temperatures relevant to collision- 
ally ionized plasmas, the earlier data were obtained using rather 
crude parameterization formulas that are not as reliable as our 
computed MCBP rate coefficients. The results reported here 
serve as what we assess to be the most accurate and complete 
set of Ar-like DR data for use in spectral diagnostic studies of 
laboratory and astrophysical plasmas. 
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Fig. 1. C omparison of existi ng total DR ground-level rate coefficients for Kn (left) and Cam (ri ght): (a) red solid curve, empirical 
results of lMewe et alj d 1980b ; (b) black solid curve, recommended data of iMazzotta et al.l d!998f) : red dashed curve, present MCBP 
results. 




kj (eV) kj (eV) 



Compa rison of existing total PR gro und-level rate coefficients for Sc iv (le ft) and Ti v (right): (a) 



Fig. 2. Compa rison of existing total 1 
experiments bv lSchippers et al. | (| l998. 



ray open circles, TSR 

2002); (b) black doted curve, compilation by Land ini & Fossil ( 19911): (c) black dash-dotted 



curve, empirical formula of Hahnl (Il99ll): (d ) black solid curve, recommended value bv IMazzotta et al.l ( 1998 ): (e) black dashed 
curve, empirical results of iMewe et al. f (fl980l) . The present RR results are also shown as the long-dashed curve. 
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Fig. 3. C omparison of existin g total DR ground-level rate coefficients for V v i (left) and Cr vn (right ): (a) red solid curve, empirical 
results of iMewe et aL ( 19801): (b ) black solid curve, recommended data of iMazzotta et alj (119981) : (c) black dash-dotted curve, 
empirical formula of Hahnl l 1991b : red dashed curve, present MCBP results. 




Fig. 4. Comparison of exist ing total DR ground-level rate coefficients. On the left is for Mn v m: (a) red solid curve, empirical 
re sults of Mewe et al.l (Il980h : (b) black solid curve, recommended data of IMazzotta et al.l (Il998h : (c) green solid curve, LS results 
of Jacobs & Davisl (119831) : red dashed curve, present M CBP results. On the right is Fe ix: (a) red solid curve, empiri cal results of 
IMewe et alJ ( 19801) : (b) blue solid curve, compil ation of[Arnaud & Raymond! (II 992h : (c) black solid curve, result of |Jacobs~et 
(1 19771) 
(0991 



(d) green solid curve, empirical data of iKato & Asanol d 19991) : (e) black dash-dotted curve, empirical formula of l Hah 
(f) "4CF" and "5CF" MCBP results of iBadnelll d2006l) : red dashed curve, present "18CF" MCBP results. The present RR 



[aD 
ahnl 



results are also shown as the long-dashed curve. 
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,,.10+ T.T-^ + 

e + Ni — > Nt 




Fi g. 5. Comparison o f existing total DR ground-level rate coeffici ents. On the left is for Co x: (a) red solid curve, empirical results 
of lMewe et al. | ( fl980l) : (b) black solid curve, recommended data of lMazzotta et al.l dl998h: re d dashed curve, present MCBP results. 
On the right is for Nixi ion : (a) red solid curve, empirical results of iMewe et alJ (11 980): (b) black solid curve, recommended 
data of Mazzotta et al. dl998l) : (c) black dash-doted curve, empirical formula of lHahrJ dl991l) : (d) blue solid curve, compilation of 
lArnaud & Rothenflugl d!985l) : (e) green solid curve, empirical data of lKato & Asanol d 19991) : red dashed curve, the present MCBP 
results. 




Fig. 6. C omparison of existin g total DR ground-level rate coefficients for C u xn (left) and Zn xm (r i ght): ( a) red solid curve, empirical 
results of IMewe et al.l dl980h : (b) black solid curve, recommended data of iMazzitelli & Mattiolil d2002l) : red dashed curve, present 
MCBP results. 
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Fig. 7. Present Maxwellian-averaged DR and RR rate coefficients for Kn (left) and Cam (right). The red and gray dashed curves 
are the total DR and RR rate coefficients, respectfully. The intra-shell An c = contributions (solid blue curve) are due to dominant 
3/7 — > 3d core excitations; the inter-shell A« c = 1 contributions (solid black curve) are due to 3p — > At' promotions. 




kT (eV) kT (eV) 



Fig. 8. Present Maxwellian-averaged DR and RR rate coefficients for Sc iv (left) and Ti v (right). The red and gray dashed curves 
are the total DR and RR rate coefficients, respectfully. The intra-shell A« c = contributions are due to dominant 3p — * 3d core 
excitations; the inter-shell A« c = 1 contributions are due to 3p — > 4-C promotions. 
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Fig. 9. Present Maxwellian-averaged DR and RR rate coefficients for V vi - Zn xm ions. The red and gray dashed curves are the total 
DR and RR rate coefficients, respectfully. The intra-shell A« c = contributions are due to 3p — > 3d (solid black curve), 3 5 — > 3d 
(solid blue curve), and 3p 2 — > 3d 2 (solid cyan curve) core excitations; the inter-shell An c = 1 contributions are due to 3p — > At' 
(solid green curve) and 3s — > AC (solid gray curve) promotions. 
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Table 1. Radial scaling parameters, A„i, for the n — 3 valen ce orbitals optimized in the Slater-type-orbital model potential. Details 
about the full optimization procedure are reported earlier bv lNikolic et al.l (12.0091) . 





K + 


Ca 2+ 


Sc 3+ - a 


-jv4+,b 


V 5+ 


Cr 6+ 


Mn 7+ 


Fe 8+ 


Co 9+ 


Ni 1,)+ 


Cu 11+ 


Zn 12+ 




1.09548 
1.11185 
1.13830 


1.09548 
1.10545 
1.16918 


1.08366 
0.97900 
1.27700 


1.05653 
0.98950 
1.17400 


1.07085 
1.14255 
0.89575 


1.07085 
1.11729 
0.87575 


1.07085 
1.09171 
0.87575 


1.05990 
1.07936 
0.84905 


1.0000 
1.0693 
0.8391 


1.0000 
1.0608 
0.8191 


1.0000 
1.0453 
0.8191 


1.0000 
1.0318 
0.8191 



adopted from: 11 Nikolic et al. (2010); b Nikolic et al. (2009). 



l±p] 



Table 2. Dominant radiative transition data from the ground state of the recombining ion 3p ( So) — > 3/5 3c/ ( P°), where v' i() 
denotes v(u) x I0 ±p . In the case of Sc 3+ and Ti 4+ ions, refer to lNikolic et alJ (l2010h and lNikolic et alJ d2009h . respectively. 



Ion 


A r 
(ns _1 ) 


S 
(a.u.) 


f 


Ref. 


Ion 


A r 
(ns _1 ) 


S 
(a.u.) 


f 


Ref. 


K + 


0.2785 [+2] 


5.004 


3.070 


Present 2 




0.836 [+2] 


5.678 


4.818 


Present 3 


n ooo/i [+2] 
U.ZZz4 L ' 


a 1 An 
4. 14 / 


Z.J 1 J 




Ca 2+ 


0.83 (15) 


^ f,A 
J.04 


A 

4. / 3(88) 






1.4610 L+2J 


2.448 


3.313 


Present 2 


0.638 L+2J 


4.614 


3.834 


MCHF C 


y5+ 


1.46 L 1 


2.45 


3.31 


ATT O Tbl* 

NIST ' 




0.682 L 1 


A 11 A 

4.210 


3.688 


CIV3 J 




1.8447 [+2] 


2.968 


4.070 


MCHF C 


Cr 


1.65^ 


2.04 


3.05 


Present 3 




l.858 L+ZJ 


1.754 


2.874 


Present 3 


l.67 L+ZJ 


2.06 


3.09 


NIST B+ 


Mn 7+ 


l.85 [+2] 


1.75 


2.86 


NIST b * 




2.05 [+2] 


2.43 


3.69 


MCHF C 




l.963 [+2] 


2.033 


3.229 


MCHF C 




2.009 L+2J 


1.489 


2.645 


Present 3 




2.190 L+ZJ 


1.300 


2.486 


Present 3 




2.0l [+2] 


1.49 


2.65 


NIST bt 


Co 9+ 


2.200 [+2] 


1.300 


2.500 


NIST b * 




2.097 [+2] 


1.724 


2.958 


MCHF c 


2.224 [+2] 


1.482 


2.726 


MCHF C 


Fe 8+ 


2.438 [+2] 


1.678 


3.054 


FAC4 d 




l.448 [+2] 


1.409 


2.628 


TDCHF f 


2.266 [+2] 


1.680 


2.983 


CIV3 e 




2.339 L+ZJ 


1.132 


2.317 


Present 3 




l.356 [+2] 


1.073 


2.843 


TDCHF f 




2.340 [+2] 


1.130 


2.320 


NIST bt 




2.889 [+2] 






SELT g 


Ni 10+ 


2.345 [+2] 


1.288 


2.528 


MCHF C 




2.404 [+2] 


1.637 


2.99 


IRON h 


l.537l +2 ] 


0.808 


2.444 


TDCHF f 




2.530 L+/J 


1.01 


2.204 


Present 3 




2.786 [+2] 
2.84l [+2] 


1.274 
1.204 


2.657 
2.575 


FAC3 h 
CIV3 1 


Cu 11+ 


2 50l [+2] 
2.462 [+2] 


1.07 (11) 
1.13 


2.29( 7) 
2.355 


NISX b,§ 
MCHF C 




2.720 L+ZJ 


0.906 


2.101 


Present 3 




l.624 [+2] 


1.08 


2.285 


TDCHF f 


Zn 12+ 


2 653 [+2] 


0.99 (11) 


2.21(7) 


NISX b,§ 












2.576 [+2] 
l.708 [+2] 


1.001 

0.958 


2.204 
2.145 


MCHF C 
TDCHF f 













assigned uncertainties: < 18%, *' < 25% 
present work: 2894-level MCBP results 
critically compil ed experimental data oflShirai et alj feOOO) 
MCHF results of llrimia &] Froese-FischeH (12003)) 
6284-level FA C results oflAggarwal et alJ ( 12006 
CIV3 results of lVerma et al. | J2006h 



f single-config. TDCHF r esults oflGhosh & Mukherjed i 19971) 

E semiempirical results of lLogino v & Tuc hkinl d2002h 

h 6164-level FAC results o flAggarwal & Keenanl d2008h 

1 restricted CIV3 results of IVerma et all d2007l) 

j 18-config. CIV3 LS results of lBaluial dl986T) 

§ extrapolated along the isoelectronic sequence 
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Table 3. The lowest An c = 0, 1 core excitation thresholds (in Rydbergs) for K n and Ca m. For brevity, we show only the core 
excitations that are below the dominant 3p 6 ( 1 Sq ) — > 3p 5 3d ( l P°) thresh old. The exception is the 3p 5 4d ( P°) t hreshold in K n 
whose configuration is marked as questionable bv lSugar & Corliss! ( 1 1 985b and only recently properly identified bv lPettersen et al.l 
(120071) . 







K + 








Ca 2+ 






K Config. 


Level(mix) 


Present 2 


NIST b 


MCHF C K Config. 


Level(mix) 


Present 3 


NIST d 


MCHF C 



9 
10 
11 

12 

13 

14 
15 
16 

17 

18 
19 

20 
21 

22 
23 
24 
25 
26 
27 



3s 2 3p 6 

3s 2 3p 5 4s 
3s 2 3p 5 4s 
3s 2 3p 5 3d 

3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3pHs 

3s 2 3pHs 

3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 

3s 2 3p 5 3d 

3s 2 3p 5 3d 

3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 

3s 2 3pHp 

3s 2 3pHp 
3s 2 3p 5 4p 
3s 2 3p 5 4p 
3s 2 3p 5 4p 
3s 2 3p s 4p 
3s 2 3p 5 4p 
3s 2 3p 5 4p 
3s 2 3p 5 4p 
3s 2 3p 5 4p 
3s 2 3p 5 3d 
3r3p 5 4d 



'S (96.3%) 
3 P° (93.2%) 

3 P° (59.5%) 
3 i^(81.1%) 
3 P° (82.4%) 

l P\ (74.8%) 

3 F° (82.4%) 
3 F 3 °(81.1%) 
3 F°(81.3%) 

3 D° (35.7%) 

'D°(41.6%) 

3 D°(71.5%) 
3 Z)°(41.9%) 
'F 3 °(34.1%) 

3 S i (97.7%) 

3 D 3 (99.9%) 
3 D 2 (75.5%) 
3 D! (76.3%) 
'02(51.0%) 
'Pi (52.5%) 
3 P 2 (55.9%) 
3 P„ (99.5%) 
3 Pi (65.5%) 
'5 (95.2%) 
'P° (90.5%)* 
'P° (90.0%) 



0.000000 0.000000 0.000000 1 



1.480381 
1.487826 
1.490116 

1.499706 
1.505571 
1.506735 

1.517409 

1.547896 
1.555741 
1.562998 

1.633040 

1.639503 

1.644348 
1.645874 
1.645999 



1.693675 
1.697258 
1.705934 
1.712827 
1.724516 
1.727938 
1.730260 
1.733330 
1.829238 
1.840715 
2.033261 



1.480834 
1.487481 
1.489303 

1.498961 
1.502934 
1.504916 



1.400644 
1.408491 
1.431384 

1.430119 
1.436726 
1.414847 



1.516872 1.437439 



1.549602 
1.556728 
1.563029 

1.634775 

1.634350 

1.644371 
1.643859 
1.645919 



1.480702 
1.488382 
1.494070 

1.561671 

1.560836 

1.570061 
1.569134 
1.571404 



1.698454 
1.701166 
1.708872 
1.714552 
1.724477 
1.728285 
1.729293 
1.732598 
1.774895 
1.840373 
2.033256 



1.618318 
1.621064 
1.628229 
1.633790 
1.643195 
1.646621 
1.648054 
1.650821 
1.702454 
1.817999 



9 
10 
11 

12 

13 

14 
15 
16 



1.663215 1.669479 1.591351 17 



3s 2 3p 6 

3s 2 3p 5 3d 
3s 2 3p s 3d 
3s 2 3p s 3d 

3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 

3s 2 3p 5 3d 

3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 

3s 2 3p 5 3d 

3s 2 3p 5 4s 

3s 2 3p s 4s 
3s 2 3p 5 4s 
3s 2 3p s 4s 

3s 2 3p 5 3d 



l S o (99.0%) 
3 P° (95.7%) 

3 F°(95.1%) 
3 F° (93.9%) 
3 F° 2 (94.0%) 

'D° (58.9%) 

3 D° (53.7%) 
3 D° (92.9%) 
3 D° (54.0%) 

'F° (57.8%) 

3 P° 2 (99.4%) 

3 P° (88.0%) 

'PI (87. 6%) 
'P?(63.1%)* 



0.000000 0.000000 0.000000 



1.879530 
1.884279 
1.893907 

1.952880 
1.963210 
1.972601 

2.086507 

2.080602 
2.093231 
2.095232 

2.081461 

2.313063 

2.327118 
2.341390 
2.369020 



1.853273 
1.857636 
1.866664 

1.934711 
1.944456 
1.953156 

2.057880 

2.062503 
2.072514 
2.072117 

2.097941 

2.210253 

2.222858 
2.238180 
2.257176 



1.760169 
1.764514 
1.773134 

1.844004 
1.853118 
1.861216 

1.972081 

1.976286 
1.985320 
1.984866 

1.995780 

2.114045 

2.126335 
2.140287 
2.160581 



2.545658 2.545658 2.492690 



present work: 51-level (K + ) and 69-level (C a 2+ ) MCBP resu lts 

UV spark spectr oscopy experimental data oflPettersen et all d2007h 
MCHF results of llrimia & Froese-Fischerl d2003 | ) 
critically compiled experimental data of lSugar & Cor liss ( 1985) 
dominant excitation threshold - see Table [2] 



Table 4. The lowest A« c 
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V 5+ Cr' 



K 


Config. 


Level(mix) 


Present 11 


NTST 


MCHF 


Level(mix) 


Present" 


NIST 


~\ A /~~~* 1 1 1 'C. 

MCHF 


1 


3s 2 3p 6 


'50(96.4%) 


0.000000 


0.000000 


0.000000 


[ S (96.6%) 


0.000000 


0.000000 


0.000000 


2 
3 
4 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 


3 P° (97.0%) 
3 P° (96.9%) 
3 P° 2 (96.7%) 


2.739060 
2.748293 
2.767020 


2.808068 
2.819413 
2.842952 


2.806866 
2.818090 
2.841350 


3 P° (97.1%) 
3 P° (96.9%) 
P° (96.6%) 


3.059906 
3.072680 
3.098672 


3.109055 
3.123582 
3. 154233 


3.110232 
3.124881 
3.155270 


5 
6 
7 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 


S F\ (97.3%) 
3 F 3 ° (96.2%) 
3 F\ (96.6%) 


2.896923 
2.906694 
2.917192 


2.941329 
2.961235 
2.981802 


f\ i~\ A A A 

2.944497 
2.963277 
2.982698 


ipo (974%) 

3 F° (95.7%) 
3 F° (96.2%) 


3.234927 
3.247762 
3.262960 


i~i f> C Ci 1 i~l Ci 

3.258178 
3.282128 
3.308454 


3.264649 
3.287197 
3.312214 


8 


3s 2 3p 5 3d 


'Z>°(71.5%) 


3.110770 


3.145142 


3.152963 


'£>° (67.4%) 


3.470531 


3.487258 


3.497713 


10 

n 


3s 3p 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 


3r\0 /TO '1/77 \ 

D" (78.2%) 
3 D<; (96.8%) 
3 £>°(71.6%) 


3.119456 
3.126389 
3.132258 


3.148578 
3.170297 
3.174174 


3.156780 
3.176499 
3.180891 


Dj (72.3%) 
3 L>° (96.8%) 
3 L>° (71.8%) 


3.478260 
3.489568 
3.499112 


3.487760 
3.515926 
3.523110 


3.499169 
3.525054 
3.532855 


12 


3s 2 3p 5 3d 


'F 3 °(77.6%) 


3.145710 


3.195308 


3.203371 


'F°(71.3%) 


3.512852 


3.546890 


3.557753 


13 


3s 2 3p 5 3d 


'/*■ (95.4%)* 


4.059132 


4.059108 


4.114305 


(95.6%)* 


4.492846 


4.492869 


4.553079 


14 
15 
16 


3s3p 6 3d 
3s3p 6 3d 
3s3p 6 3d 


3 £>i (69.1%) 
3 D 2 (69.1%) 
3 D 3 (69.2%) 


5.013754 
5.017989 
5.024692 


5.007759 
5.010726 
5.015474 


5.122273 
5.124985 
5.129131 


3 D 1 (71.9%) 
3 D 2 (71.8%) 
3 Z) 3 (72.0%) 


5.547059 
5.552838 
5.562141 


5.546697 
5.550917 
5.557707 


5.667028 
5.670942 
5.676969 


17 


3s3p 6 3d 


l D 2 (63.6%) 


5.191350 


5.161717 


5.318685 


t> 2 (66.4%) 


5.742414 


5.721178 


5.883764 



a present work: 2894-level MCBP results 

b critically compil ed experimental data oflShirai et all d2000h 

MCHF results of llrimia & Froese-F ischer (20Q3|) 

* dominant excitation threshold - see Table [2] 
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Table 5. The lowest An c = 


core excitation thresholds (in 


Rydbergs) for Mn vm and Fe ix. 
















Mn 7+ 












Fe 8+ 








K 


Coring. 


Level(mix) 


Present" 


NIST b 


MCHF 


Level(mix) 


Present" 


NISI* 


MCHF C 


CIV3 e 


IRON 1 


FAC3 g 


l 


3.r3p 6 


'So (96.7%) 


0.00000 


0.00000 


0.00000 


'5(1(96.9%) 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.0000 


2 
3 
4 


3s 2 3p s 3d 
3s 2 3p !i 3d 
3s 2 3p 5 3d 


3 P° (97.2%) 
3 P° (96.9%) 
3 P° (96.4%) 


3.36224 
3.37921 
3.41384 


3.40502 
3.42372 
3.46275 


3.19510 
3.21381 
3.25264 


3 P° (97.4%) 
3 P° (96.9%) 
3 P° (96.2%) 


3.67368 
3.69563 
3.74058 


3.69767 
3.72084 
3.76963 


3.47859 
3.52725 


j.Oy 1 DO 

3.72078 
3.76962 


J. /OJ04 

3.78883 
3.84039 


j. IH 1 J 

3.7710 
3.8200 


5 
6 
7 


3s 2 3p 5 3d 
3s z 3p 5 3d 
3s 2 3p 5 3d 


3 F° (97.5%) 

ini i c\ a c\rn \ 

{•" (94.9%) 
3 F\ (95.5%) 


3.55446 
3.571 13 
3.59269 


3.57067 
3.59878 
3.63198 


3.36662 
3.39298 
3.42463 


3 F° (97.6%) 

'1m /m r\r~i7 \ 

r" (93.9%) 
3 F\ (94.6%) 


3.88447 
3.90343 
3.93169 


3.88026 
3.91217 
3.95325 


3.64440 
3.67440 
3.71373 


3.88027 
3.91217 
3.95329 


3.97113 
4.00201 
4.04308 


3.9468 
3.9797 
4.0215 


8 


3s 2 3p 5 3d 


3 D° (68.7%) 


3.81841 


3.82162 


3.62227 


3 D° (64.5%) 


4.16736 


4.15184 


3.92030 


4.15171 


4.25326 


4.2373 


9 


3s 2 3p 5 3d 


'D° (64.4%) 


3.81303 


3.82565 


3.62467 


'Z)°(62.1%) 


4.16451 


4.16224 


3.92846 


4.16220 


4.26847 


4.2540 


10 
1 1 


3s 2 3p 5 3d 
3s 2 3p s 3d 


3 D° (96.8%) 
3 D° 2 (64.7%) 


3.83602 
3.85013 


3.85773 
3.86961 


3.65585 
3.66842 


3 £>° (96.8%) 
3 D° (62.6%) 


4.19143 
4.21227 


4.19744 
4.21567 


3.96326 
3.98226 


4.19745 
4.21573 


4.29702 
4.31956 


4.2836 
4.3039 


12 


3s 2 3p 5 3d 


'F°(67.1%) 


3.86567 


3.89595 


3.69578 


l F° (62.2%) 


4.22812 


4.24494 


4.01241 


4.24509 


4.34420 


4.3281 


13 


3s 2 3p 5 3d 


'P°(95.7%)* 4.91368 


4.91368 


4.76446 


'/» (95.9%)* 


5.32660 


5.32678 


5.14621 
5.212 d 


5.32645 


5.47936 


5.4598 


14 
15 
16 


3s3p 6 3d 
3s3p 6 3d 
3s3p 6 3d 


3 D, (73.5%) 
3 D 2 (73.4%) 
3 Z) 3 (73.7%) 


6.05426 
6.06173 
6.07394 


6.08432 
6.09007 
6.09935 


5.99785 
6.00326 
6.01167 


3 Di (74.7%) 
3 D 2 (74.6%) 
3 D 3 (74.9%) 


6.61242 
6.62228 
6.63882 


6.62249 
6.63001 
6.64254 


6.50906 
6.51629 
6.52770 


6.62249 
6.63001 
6.64254 


6.75891 
6.76762 
6.78194 


6.7176 
6.7251 
6.7377 


17 


3s3p 6 3d 


>D 2 (68.2%) 


6.26561 


6.27726 


6.23324 


'D 2 (69.4%) 


6.84122 


6.83333 


6.76228 


6.83334 


6.98165 


6.9447 



present work: 2894-level MCBP results 

critically compil ed experimental data oflShirai et alj f2000). 
MCHF results of llrimia & Froese-FischeJ d2003h . 

dominant excitation threshold - see Ta ble [2] 

single configur ation TDCHF result s of lGhosh & Mukheried dl997l) . 

CIV3 results of IVerma et all j2006h. 

IRON Project results of lStorev et al.ld2002h. 
6284-level FAC results of lAggarwal et alj 00061) . 
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Table 6. The lowest An c = core excitation thresholds (in Rydbergs) for Co x and Ni xi. Uncertainties are enclosed in lower 
parentheses. 



Co 9+ Ni 



K 


Config. 


Level(mix) 


Present 21 


NIST b 


MCHF C 


Level(mix) 


Present" 


NISI* 


MCHF C 


FAC3 e 


CIV3 f 


1 


3.r3p 6 


l S o (97.0%) 


0.00000 


0.00000 


0.00000 


'50(97.1%) 


0.00000 


0.00000 


0.00000 


0.0000 


0.00000 


2 
3 
4 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 


3 P° (97.5%) 
3 P° (96.9%) 
3 P° 2 (95.8%) 


3.97623 
4.00373 
4.06022 


3 999+ 

4.026L 
4.084^, 


3.70893 
3.73778 
3.79773 


3 P° (97.6%) 
3 P° (96.8%) 
3 P° 2 (95.2%) 


4.27882 
4.31288 
4.38307 


4.27667 
4.31002 
4.38274 


3.95697 
3.99197 
4.06469 


4.2998 
4.3349 
4.4086 


4.34736 
4.38366 
4.45723 


5 
6 
7 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
3s 2 3p 5 3d 


3 F° (97.7%) 
3 F° (92.8%) 
3 F° (93.2%) 


4.20594 
4.22715 
4.26350 


4.199* 

4 235' 
'■"■'pi] 

4.283+ 25) 


3.91676 
3.94997 
3.99798 


3 F° (97.8%) 
3 F°(91.3%) 
3 F°(91.3%) 


4.52944 
4.55153 
4.59721 


4.49482 
4.53374 
4.59342 


4.18448 
4.22023 
4.27778 


4.5391 
4.5785 
4.6393 


4.64639 
4.68656 
4.73636 


8 


3s 2 3p 5 3d 


3 D° (62.3%) 


4.50743 


4.494j 38) 


4.21215 


3 D° 3 (60.3%) 


4.84888 


4.80666 


4.49868 


4.8593 


4.93579 


9 


3 t 2 3n 5 3^/ 




4 50846 


4 51 1 1 " 

J (29) 


4.22855 




4 85545 


4.83728 


4.52643 


4.8953 


4.97061 


i n 
11 


JJ JU JCl 

3s 2 3p 5 3d 


Ls y \~7yj. 1 /v ) 

3 D° (61.4%) 


4 S^Q4^ 

4.56821 


4 S^fil I 

4.573; 24) 


4 766^4 
4.29382 


} n° cofi s^m 
3 D° (60.2%) 


t.0077J 

4.92948 


4 87^7^ 
4.91337 


4 Sfifil 7 

t.JUU 1 z. 

4.60454 


4 Q7Q0 

4.9687 


S 07SQQ 

J.UiJ77 

5.04846 


12 


3s 2 3p 5 3d 


'F 3 ° (59.0%) 


4.58588 


4.606; 25) 


4.32726 


^ (56.1%) 


4.94891 


4.95018 


4.64195 


5.0024 


5.12673 


13 


3s 2 3p 5 3d 


'P°(96.1%)* 


5.73579 


5.73579 


5.51933 
5.5952 d 


(96.2%)* 


6.14044 


6.14158 


5.88618 
5.9728" 


6.2574 


6.41884 


14 
15 
16 


3s3p 6 3d 
3s3p 6 3d 
3s3p 6 3d 


3 Dt (74.2%) 
3 D 2 (74.0%) 
3 D 3 (74.4%) 


7.21441 
7.22694 
7.24842 


7.178+ 
7.186* 
7.198+ 31) 


7.01984 
7.02926 
7.04436 


3 £>,(75.1%) 
3 D 2 (74.1%) 
3 £> 3 (75.3%) 


7.76928 
7.78502 
7.81290 


7.720+ 
7.730* 
7 740 r 

'■""(33) 


7.53153 
7.54354 
7.56317 


7.7562 
7.7687 
7.7904 


7.85347 
7.86852 
7.89090 


17 


3s3p ( '3d 


'£> 2 (68.8%) 


7.46032 


7.410+33, 


7.29060 


l D 2 (69.7%) 


8.03277 


7 970 r 

'•""(34) 


7.81976 


7.9965 


9.00447 



a present work: 2894-level MCBP results 

b critically compil ed experimental data oflShirai et alj d2000h . 

c MCHF results of llrimia & Froese-F ischer (2003]). 

+ extrapolated along the isoelectronic sequence. 

* dominant excitation threshold - see Ta ble [2] 

d single-configuration TDC HF results of iGhosh & Mukherieel jl997l) . 

e 6164-level FAC results o flAggarwal & KeenarJ d2008h . 

f restricted CIV3 results of IVerma et alj d2007l) . 



16 



D. Nikolic et at Dielectronic Recombination of Argon-Like Ions 



Table 7. The lowest An c = core excitation thresholds (in Rydbergs) for Cu xn and Zn xm. Uncertainties are enclosed in lower 
parentheses. 



Cu 11+ 



Zn' 



K 


Config. 


Level(mix) 


Present' 1 


NIST 


MCHF 


Level(mix) 


Present 3 


NIST 


MCHF 


I 
i 




o \y I .z* /c ) 


v. www 


n nnnnnn 
v. www 


n nnnnnn 
v. www 


a q \y 1 .J /c ) 


n nnnnn 
v.wvw 


n nnnnn 

V.VWW 


n nnnnn 


2 
3 
4 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
5s 5p id 


3 P° (97.7%) 
3 P° (96.6%) 

3r>0 /O /I yl /97 \ 


4.56676 
4.60819 
4.DV3 /o 


4.638* 32) 
4.711 (25) 


4.19956 
4.24146 
4.32844 


3 P° (97.7%) 
3 P° (96.3%) 

3 DO fCYl A (77 \ 

r 2 (yjA /c) 


4.85185 
4.90147 
5.00410 


4 905 § 
4 943 § 


4.43697 
4.48650 
4.58924 


5 
6 

7 


3 5 3/? 3a 
3s 2 3p 5 3d 

JA J^/ Jtl 


3t7'0 i'CVH C\<7! \ 

t '4 (97.9%) 
3 F° (89.6%) 


a r> r> 

4.83899 
4.86273 

4. Q901 S 


4.837 (36) 
4.880 ( § 35) 

4 Q4.1 § 


AAA Or\H 

4.4480 / 
4.48537 
4 ss^nq 


3r70 /m c\(77 \ 

i*4 (97.9%) 
3 F° (87.7%) 

3 F° CR4 Q%\ 


5.14/10 
5.17130 

S 94.1 Sf> 


c 1 c^.§ 

j.ijo (38) 
5.203 (38) 
s ?7n^ 

3.Z/U (31) 


4./0801 
4.74554 
4 89^sn 


8 


3s 2 3p 5 3d 


! D° (57.2%) 


5.19405 


5 203 ! 


4.82348 


l D° (55.4%) 


5.53462 


5.548f 36) 


5.12106 


9 
10 
1 1 


3s 2 3p 5 3d 
3s 2 3p 5 3d 
3i 2 3n 5 3d 


3 D° (59.8%) 
3 Z)° (96.2%) 
3 D° f 59 4%") 


5.17774 
5.23130 
5 28354 


5 177 § 

'(44) 

5.21500 

5 280 § 


4.78044 
4.86330 
4.91565 


3 Z)° (59.6%) 
3 £)° (95.8%) 
3 D° (58 7%") 


5.50490 
5.57297 
5.64097 


5 519 § 
5.55700 
5 633 § 


5.05792 
5.15855 
5.22837 


12 


3s 2 3p 5 3d 


'F° (54.3%) 


5.30733 


5.320 (29) 


4.95778 


l F° (52.7%) 


5.67000 


5.676 ( § 31) 


5.27605 


13 


3s 2 3p 5 3d 


! P° (96.2%)* 


6.54764 


6.54764 


6.24830 
6.346 d 


'P° (96.1%)* 


6.95380 


6.95379 


6.60700 
6.7154 d 


14 
15 
16 


3s3p 6 3d 
3s3p 6 3d 
3s3p 6 3d 


3 £>! (75.6%) 
3 £> 2 (74.9%) 
3 £> 3 (75.9%) 


8.29610 
8.31511 
8.34975 


8 - 268 (27) 

8.277f 36) 
8.29lf 36) 


8.04518 
8.06020 
8.08529 


3 £>! (75.9%) 
3 D 2 (75.1%) 
3 £> 3 (76.3%) 


8.82751 
8.85018 
8.89283 


8.814j 29) 
8.822 ( s 38) 
8.838 ( § 39) 


8.56194 
8.58042 
8.61207 


17 


3s3p 6 3d 


^2 (70.1%) 


8.57720 


8.538 ( § 39) 


8.35116 


'£> 2 (70.3%) 


9.12726 


9 102 § 


8.88624 



a present work: 2894-level MCBP results 

b critically compil ed experimental data oflShirai et alj d2000T ) 

MCHF results of llrimia & Froese-F ischer ( 20Q3|) 

s extrapolated along the isoelectronic sequence 

* dominant excitation threshold - see Ta ble [2] 



single configuration TDCHF results of Ghosh & Muk heried i 19971) 
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Table 8. Fitting parameters E, (K) and c, (cm 3 s 'K 3 ^ 2 ) used for modeling the Maxwellian-averaged DR rate coefficients - see 
Eq. (0. Uncertainties are enclosed in lower parentheses, where v^F denotes v(u) x \Q ±P . 



ion 




Ei 


Ei 


E } 


Es 


Et 


Ej 


Ci 


Cl 


C3 


c 4 


C5 


C6 




K + 


2.451™ 


3.504™ 


4.094™ 


4.766™ 








6.292™ 


2.350™ 


1.165™ 


2.470™ 








Cir 


2.282™ 


3.682™ 


4.479™ 










3.843™ 


8.040™ 


8.670™ 










Sc 3+ 


Lies;; 5 ' 


2.188™ 


4.771™ 


6.814™ 








9.817™ 


6.241™ 


2.459"' 


6.990j~ 3] 








Ti 4+ 


1.5603JJ, 11 


9.8857™ 


3.7145™ 


1.4089™ 


5.9207™ 


4.1099™ 


5.1296™ 


1.006™ 


1.2449[7 ( 61 


4.3180™ 


1.8210™ 


7.9042™ 


6.040™ 


8.457™ 


V 5+ 


1.655,'-' 


4.414™ 


1.088™ 


2.299™ 


4.076™ 


6.415™ 


7.743™ 


1.968™ 


1.996™ 


3.053™ 


2.020™ 


1.877™ 


1.031™ 


1.569™ 


Cr 6+ 


5.649,™ 


1.664,™ 


5 521 1+41 


1.263™ 


2.463,™ 


5.499™ 


7.764™ 


3.022™ 


2.606™ 


4.024™ 


2.740™ 


1.102™ 


5.602™ 


1.199™ 


Ma 1 * 


3.284™ 


Ull™ 


2.818] 4 * 1 


7.017™ 


1.681™ 


4.343™ 


8.056™ 


4.596™ 


1.499™ 


1.960™ 


5.060™ 


7.040™ 


3 42l'- 2J 


1.963™ 


Fe s+ 


2.02 1™ 


2.187™ 


7.063™ 


2.074™ 


4.974™ 


2.743™ 


8.348™ 


3.786™ 


1.494™ 


9 447™ 
'■"'in 


2.100™ 


3.520™ 


2.785™ 


2.734™ 


Co ,+ 


1.277™ 


9.012™ 


2.895,™ 


9.202™ 


1.782™ 


4.485™ 


9.033|2! 51 


1 549™ 


2.685™ 


4.341™ 


3.670J;, 31 


1.555™ 


5.304™ 


2.857I;, 11 


Ni 10+ 


1.929™ 


9.502™ 


2.690™ 


6.006,™ 


1.438™ 


4.627™ 


1.020™ 


1.715™ 


5.967™ 


2.660™ 


7.130™ 


1.782™ 


8.587™ 


3.393™ 


Cu" + 


4.706™ 


1.541™ 


2.982™ 


7.313™ 


2.166™ 


6.000™ 


1.160™ 


5.843™ 


1.690™ 


6.440™ 


1.388™ 


2.954,™ 


1.494™ 


3.518™ 


Zn 12+ 


1.182™ 


5.335™ 


1.502™ 


4.288™ 


1.626™ 


4.313[+ 51 


1.034™ 


9.08157^ 


1.299™ 


1.492™ 


3.050™ 


3.209™ 


1.049,™ 


4.396™ 
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Table 9. Fit coefficients for total ground state RR rate coefficients of recombining ions, see Eq. ©, where vf , denotes v(u) x 10'' 
with uncertainties given in lower parentheses. 



ion a 


A 

(1CT 11 cm^" 1 ) 


B 


To 
(K) 


T\ 
(K) 


L 


T 2 
(K) 


K + 


4.528[+ 0] 


4 234 [_1] 


5.931™ 


2.897™ 


3.049™ 


1.645™ 


Ca 2+ 


2.248[ 4 + } 1] 


6.605™ 


6 175 [+0] 


6.032;; } 6] 


3.15 8^2^ ^ 


2.100™ 


Sc 3+ 


4.196™ 


2.250[" 1] 


4.232[+ 2] 


1 110 [+8] 


3.993™ 


3.973™ 


Ti 4+ 


3.989™ 


5.658[- 1] 


9 517 [+1] 

y.jL / (3) 


4.776™ 


1.065™ 


3.534™ 


y5+ 


4.193™ 


5.887[~ 1] 


1 993 [+2] 


2.207™ 


9.247™ 


2.012™ 


Cr 6+ 


5.744™ 


6.220[~ 1] 


2.282[+ 2] 


2.25o; 8 + ; ] 


4.925™ 


2.948™ 


Mn 7+ 


2.976™ 


4.838™ 


1 397 [+3] 


4.566™ 


5.863™ 


1 672 [+4] 

i.u/z- (3) 


Fe 8+ 


3.341™ 


4.865|~ 1] 


1.891™ 


5.m™ 


5.747™ 


2 734 [+4] 


Co 9+ 


2 542 [+1] 


3.600}" 1] 

(o) 


5.026™ 


5.480™ 


1 142™ 


2.509™ 


Ni 10+ 


Till 

3.076; 2 + ; ] 


r in 

3.849[~ 1] 


5.284|; 3] 


4.945J* 71 


r 11 

1 021 L J 


r ■ a i 

2.571™ 


Cu 11+ 


3.527J+ 11 


4.052; 3 - 1] 


5.855^ 3 ^ 


5.691™ 


8.842™ 


3.626™ 


Zn 12+ 


3.964™ 


4.107| 3 - 1] 


6.576[+ 3] 


5.855^ ^ 


8.595™ 


3 942 [+4] 



present work: t < 200, 3 < n < 1000 



